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1. Introduction 

Lattice QCD has widely been used to get information about the properties of matter at high 
temperature and vanishing net baryon density [jl]]. These studies have been extended to non- 
vanishing baryon density. I.e., the equation of state has been discussed using Taylor expansion [§], 
reweighting ||3|] or imaginary chemical potential techniques 

Furthermore the free energy of static quark anti-quark sources at ju& = received much atten- 
tion |@. We will now apply here the Taylor expansion approach in order to extend the results on 
this free energy to non-zero quark chemical potential. A first attempt to do this was discussed in 
[^]. Furthermore we will analyze the large distance behaviour of the free energy and determine the 
corresponding expansion coefficients of the screening mass, which can be compared to perturbative 
predictions. 

We will restrict the discussion of the expansion up to the 2nd order in here. Results up to 
the 6th order as well as furhter details on the singulation and the Taylor expansion technique can 
be found in 

2. Taylor expansion method in QCD at finite density 

The staggered fermion partition function 

= JDU A(h) e- s , (2.1) 

where in our 2-flavour case A(ju) is the square root of the fermion determinant, can be Taylor 
expanded in powers of the quark chemical potential /I = jU&/3 by inserting the expansion of A(ju) 
in powers of jl. 

A(/i) = A(0) (1 +D { pL+D 2 pL 2 + •••). (2.2) 

The odd orders in the expansion of Z^ are vanishing because D n is imaginary for odd and real for 
even n. The simulation can be done at pL = 0, where the D n are handled like observable quantities. 
This procedure can be extended easily to the /J. -dependent expectation value of an observable 
that does not directly depend on /x, 

K), = — [dU SAM e- S = ^•l)o+(^)oM + -- 4 +(^y 7 

= {£/) Q (l+a l n+a 2 H 2 + --- + a 6 n 6 ) . (2.3) 

We apply this scheme to the Polyakov loop correlation functions discussed in the next section. 

3. Heavy quark free energies 

A heavy (static) quark Q at site x is represented by the Polyakov loop, 

L{x) = Y[ U 4 (x,x 4 ) , (3.1) 

X4 = \ 
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Figure 1: The order coefficients Q and Q for the singlet (a) and colour averaged (b) free energies. 
Q is matched to the T = heavy quark potential at small distances (a). 

which is an SU(3) matrix. A heavy anti-quark Q is described by the corresponding hermitian 
conjugate matrix. The colour averaged and singlet QQ-free energies F^^r, T,n) and Fgg(r, J \fi) 
can be calculated from the corresponding correlation functions <^ av,1 (r), 

F^\r,T,n) = -T \n{^\r))=f^(rj)+f^ 2 (rj) + & ((£) 4 ) , (3.2) 

where 

^ y ™c x,y 
■ Jy iy c x,y 

In Fig. [I] and || we show the results for the Oth and 2nd order expansion coefficients in fi/T for 
temperatures above and below T c . Because f^A 2 ( r ^ T) is always negative, we find that for small /I 
the free energy of a static quark anti-quark pair decreases in a medium with a net excess of quarks 
or anti-quarks. 

We determine the large distance value of these coefficients, shown in Fig. |3L by taking the 
weighted average of the values at the five largest distances. f^A (°°,T) signs the transition tem- 
perature by showing a pronounced peak at T c . The coefficients of the colour averaged and singlet 
free energy are identical at infinite distance. 

4. Screening masses 

For temperatures above T c and large distances r the heavy quark free energies are expected to 
be screened, 

AF*f(r,T,n) = /S 1 (- ) r i / I )-*S 1 (r,r,M)~^^ (4.1) 
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Figure 2: The 2 order coefficients of the singlet (a,c) and colour averaged (b,d) free energies below (a,b) 
and above (c,d) T c . 



with n = 1 , 2 for the singlet and colour averaged free energies respectively. In the infinite distance 
limit we thus get, 



in 



av.l 



(T,^) = -hm- r ln( K AF^ l (r,T,^ 



(4.2) 



Expanding the logarithm in powers of pL we see that only the even orders of the expansion coeffi- 
cients for the screening mass are non-zero, i.e. the second order coefficient can be written as 



m 2 (T) 



hm , 

r-w> r Af v - ( r T) 



(4.3) 



where we write ^fgQ n for the expansion coefficients of AFqq ■ We use this expression to determine 
the expansion coefficient of order 2 from the infinite distance limit. For the 0th order coefficient we 
refer to fits of the form of eq. 1.1. In fact the rational expression under the limit in eq. 4.3 shows 
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Figure 3: The coefficients for the singlet and colour averaged free energies at infinite distance rT versus 
temperature. 



only little r-dependence and we see a wide plateau starting close to rT = 0. We therefore determine 



m,' av (T) by fitting the ratio in eq. 4.3 to a constant in the plateau range. The resulting coefficient, 



shown in Fig. |] (b) is always positive and leads to an enhancment of the screening mass for small 
/ifo ^ 0. At high temperatures m^ m '(T) is in remarkable agreement with the LO high temperature 
perturbative prediction resulting from the pL -expansion of the Debye mass 



where we introduce an additional scale factor A via m l {T) = A ■ m D (T), which is determined from 
fitting ml{T). We find A = 1.397(18). Finally we see, that m\(T) and mf(T) differ by a factor 
1 /2, which is expected at high temperature for all the coefficients but has not been observed for the 
0th order so far. 



5. Conclusions 

We analyzed the dependence of heavy quark free energies and screening on the baryon- 
chemical potential. The free energies are decreasing to first order in while the screening mass 
is increasing. This suggests that the screening length in a baryon or anti-baryon rich quark gluon 
plasma decreases with increasing value of the chemical potential, which is consistent with the ex- 
pectation that a non-zero \i shifts the transition to lower temperatures. The screening behaviour 
is in good agreement with perturbation theory for T/T c >2. We observed that the /i-dependent 
corrections of the colour averaged screening mass are twice as large as those of the colour singlet. 
This suggests that the contribution to the /I -dependent corrections of the colour averaged screening 
mass is due to two-gluon exchange. 
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Figure 4: Oth and 2nd order expansion coefficients of the screening mass. The solid lines are connected to 
the LO perturbation theory. 
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